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Abstract 



The Bonnor-Swaminarayan solutions are boost-rotation symmetric space-times which de- 
scribe the motion of pairs of accelerating particles which are possibly connected to strings 
(struts). In an explicit and unified form we present a generalised class of such solutions with 
a few new observations. We then investigate the possible limits in which the accelerations 
\Q • become unbounded. The resulting space-times represent spherical impulsive gravitational 

waves with snapping or expanding cosmic strings. We also obtain an exact solution for a 
snapping string of finite length. 
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fc-* ■ 1 Introduction 

We have recently shown that expanding impulsive spherical gravitational waves may be 
considered as impulsive limits of the Robinson-Trautman family of vacuum type N space-times. 
Such waves may be explicitly constructed using Penrose's method of cutting Minkowski (or 
de Sitter or anti-de Sitter) space along a null cone and re-attaching the two parts with a suitable 
warp. However, particular solutions with these properties may also be obtained as limiting 
cases of certain known solutions which have boost-rotation symmetry. (The first example of this 
procedure was presented by Bicak and Schmidt at the end of section 5 of j|.) The purpose of 
the present work is to list some particular boost-rotation symmetric solutions in a suitable form, 
and to investigate their possible impulsive wave limits in detail. 

Boost-rotation symmetric space-times are the only explicitly known exact solutions of Ein- 
stein's vacuum field equations which describing moving particles, are radiative and asymptoti- 
cally flat in the sense that they admit global, though not complete, smooth null infinity, as well 
as smooth spacelike and timelike infinities. The general properties of such solutions have been 
extensively investigated and reviewed by Bicak and Schmidt § (see also [|| and ||), and need 
not be repeated here. 

Here, we will concentrate on the Bonnor-Swaminarayan (BS) solutions |7]] which are specific 
boost-rotation symmetric solutions. These are particularly appropriate to consider in the null 
limit in which impulsive waves will arise. However, some related solutions also have similar 
limits. We will therefore review the BS and related solutions in a unified way in the next two 
sections, adapting the notation to our requirements and making a few new observations. We 
will then investigate their impulsive limits in section 4. 
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2 The Bonnor— Swaminarayan solutions 



The BS solutions M can be described by the line element 



ds 2 



where 
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For this line element, the metric functions ti and A depend only on p 2 and £ 2 — r 
dependence exhibits explicitly the rotation and boost symmetry. 

The detailed physical interpretation of this space-time was described in j7j and [|| (a summary 
is given in In general, it represents the motion of two pairs of uniformly accelerating particles 
which are possibly connected to conical singularities (strings or struts). The particles are located 
symmetrically, two on the positive and two on the negative £-axis. However, the point "masses" 
are not black holes, but are of the type that are usually described as "Curzon-Chazy" particles. 
The radiative properties of the BS solutions have been extensively studied by Bicak |K| |ll|, ||, 0]. 




Figure 1: A space-time diagram of the BS solutions 

The above solution contains five arbitrary constants mi, m 2 , A%, A 2 and B. These determine 
the mass and acceleration of each particle and the singularity structure on the axis of symmetry 
p = 0. Each pair contains two particles with masses rrii. These are located on the axis at the 
points where Ri = 0. They thus have hyperbolic worldlines in the t-Q plane given by 



in, 



p = o, c 




(4) 



Each particle is uniformly accelerated with respect to a Minkowski background, the accelerations 
being ±A» in the ±£ directions (we shall assume A 2 > A\ > 0). The minimum distance of 
each particle from £ = is given by . 

The metric (||) will generally contain conical singularities, representing strings or struts, 
located on the axis of symmetry. These will occur unless A + /i— >0asp— >0 (see [||]). In fact, it 
is possible to choose the constants mi, Aj and B such that this regularity condition is satisfied 
on appropriate sections of the axis. 

The following cases are of particular interest: 



Case 1 



It is possible to make the choices 



mi 



\A\A\ ' m2 \A\A\ ' 

With this, the axis is regular everywhere except at the locations of the point particles. 
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Figure 2: 



mi 


m 2 


m 2 


mi 


+ 






+ 










-Ai 


-A 2 


A 2 


Ai 



c 



This case had previously been described (in the static region) by Bondi []T^] in 1957. For this 
choice of constants, the outer particles have positive mass, and the inner particles negative mass. 
The interactions between each pair of particles causes them both to accelerate towards infinity. 

Case 2 

The conditions for Case 1 can be relaxed by taking 

(A[ 
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A\f 



AA\A\ 



m 2 arbitrary, B = 0. 
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In this case, the axis is only regular between the particles of each pair and between the two pairs 
of particles. However, the outer particles are connected to infinity by two semi-infinite strings. 



Figure 3: 
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In this case, the outer particles must have positive mass and are "pulled" towards infinity by the 
tension in the string. However, the inner particles have arbitrary mass which may be positive, 
negative or zero. 

Case 2a 

Of particular interest is the special case of 2 in which m 2 = 0. For this choice, the inner 
particles disappear, and the two remaining particles are pulled towards infinity by semi-infinite 



strings. This has been described by Bicak, Hoenselaers and Schmidt [13]. 

P 

mi 
± 

Figure 4: 
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In this particular subcase, however, the restriction that the outer particles must have positive 
mass no longer occurs. We only require 



mi arbitrary, 



m 2 = 0, 



B = 0, 



(7) 



so that the remaining particles may also have negative mass, in which case the conical singularity 
corresponds to a strut rather than a string. (It may be noticed that this subcase can also be 
obtained in the limit A 2 — * A\, in which case mi — > and the masses at the end of the conical 
singularity are replaced by the original inner masses m 2 which may be arbitrary.) 



Case 3 



A further case can be obtained by imposing the conditions that the axis is regular only 
between each pair of particles and between each outer particle and infinity. In this case there 
remains a string between the two inner particles which must have negative mass. However, the 
outer particles may have arbitrary (positive or negative) mass. The parameters are given by 
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Figure 5: 
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Case 3a 

The special case of 3 arises if mi = 0. This would result in two particles of negative mass 
connected by a string. However, the restriction on m 2 no longer occurs and it may be chosen 
arbitrarily (as in case 2a above). This has been described in [13|. Alternatively, we may obtain 
the same situation by taking the limit as A 2 — > A± in (0). This subcase is thus given by 



mi arbitrary, 



m 2 = 0, 
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Figure 6: 
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Notice that for the special cases 2a and 3a describing only two accelerated particles, the 
metric functions (|2|) can be written in a simple form (see p6|) 
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In these expressions, the unnecessary indices have been omitted (i.e. m = mi, A = Ax,R = Rx)- 
Case 4 



A further case has been described in detail in [14| for which the axis is regular everywhere 
except between the two pairs of particles. The parameters are here given by 

Ax 



rax arbitrary, 
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In this case the masses of the two particles in each pair must be of opposite sign. If the outer 
particle has positive (negative) mass, the singularity corresponds to a strut (string). 



Figure 7: 
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3 The Bicak— Hoenselaers— Schmidt solutions 



Let us now consider some further solutions that can be obtained from those above by taking 
suitable limits. A further solution will also be included which contains the special cases above. 



Case 2b 

One special limit (described in detail by Bicak, Hoenselaers and Schmidt fi~3|| ) of case 2 above 
that is of particular interest is that in which A± — ► 0. In this case, the outer particles, and the 
strings attached to them, are scaled out to infinity. The significance of this solution is that the 
remaining particles move freely and are not connected to strings. However, they move under 
the action of an exterior gravitational field for which the source is at infinity. It may be noted 
that although mi ~ A^ 3 in this limit, Ri ~ A^ 2 and p, is bounded. The resulting solution is 
given by 
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Figure 8: 
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Case 4a 

We can also consider the limit of case 4 in which the A 2 — ► A\ so that the two particles in 
each pair approach each other. This procedure has been described by Bicak, Hoenselaers and 
It is convenient to introduce the parameter Mqi 



Schmidt g4|. 
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obtained unless the parameter mi is rescaled in such a way that Mqi is kept constant, 
particular solution is given by 
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It may be noticed that an analogous limit of case 1 does not exist as there is no freedom to 
rescale the mass m± (see (||)) and Minkowski space is obtained as the two particles coalesce. 



Figure 9: 
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As shown in the metric describes a combination of monopole and dipole terms, so that 
the resulting particle has been referred to as a Curzon-Chazy (Ol)-pole particle (for this reason 
we have denoted the parameter by Mqi). The space-time simply contains two such particles. 
The situation looks similar to that for case 2b above. However, in this case there is no external 
fields and the particles are accelerated by their internal dipole component. 



Multipole particles 

In [Q other boost-rotationally symmetric spacetimes generalising the Bonnor-Swaminarayan 
solution (jl||2|) were found. These represent the fields of two accelerating particles with arbitrary 
multipole structure attached to conical singularities as in the cases 2a and 3a above. Moreover, 
for a special choice of parameters the space-times may be free of conical singularities as in the 
case 4a, or 2b but without an external field. 

This class of solutions was presented in [14| using prolate spheroidal coordinates related to 



the Weyl form of the boost-rotational symmetric metric. We will give here an explicit closed 
form of these solutions in terms of the metric functions for the line element (p]) . However, the 
expression for A is presented here in a much simpler form than that given in [Q . The derivation 
of this new expression, and the relation between the two forms, are contained in the appendix. 
In addition, using the scaling property of the Weyl metric, we can include the acceleration 
parameter A explicitly (this parameter does not appear in fl4|| , in which the scaling was used 
to put A = g)- In this new form, the general class of multipole solutions is given by 
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where the constants M n represent the multipole moments, the argument of the Legendre poly- 
nomials P n is a = (1 — xy)/(x — y), and C, D are constants to be specified below. For the 
prolate spheroidal coordinates x and y and for a in (|i~5|) one has to substitute from the relations 



x-y = AA 2 R, x + y = 2A 2 (p 2 + ( 2 -T 2 ), 

where R is given by (jll]). These imply the useful formulae 

(x 2 - 1)(1 - y 2 ) = 16.4V (C 2 - r 2 ), 1 + xy = 2A 2 (t 2 - r 2 - p 2 ). 

It now remains to set the values of the constants C and D. To ensure that the space-time is 
regular on the "roof" ( 2 — r 2 = 0, it is necessary that p = A at p = = ( 2 — t 2 (see 
This yields explicitly that D = C + Y^n=o M n /2 n . Moreover, the metric is regular on the axis 
provided p + A = at p = 0. Regularity of the axis between the two particles, analogous to 
the case 2a (see Fig. 4), requires that 

00 M 

C = -£^, D = 0. (17) 

n=0 Z 

In this case there is generally a string connecting the particles to infinity. The alternative 
situation, analogous to case 3a (see Fig. 6), in which there is a string between the particles and 
the axis is regular outside is given by 

C = 0, D = J2^. (18) 

n=0 Z 

The axis is obviously regular everywhere, except at the particle, if the combination of multipole 
moments satisfies the condition that J2^ =0 M n /2 n = so that C = = D. This was pointed 
out in 



Considering only the case n = in (15), we recover the previous formula (|To| ) for the monopole 
BS particles as in case 2a or 3a with identification Mq = —4mA. 

It is also straightforward to write an explicit solution representing accelerated dipole particles 
with the single moment Mi'. 
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With these two observations, it may also be seen that case 4a given by (|l4|) is a special case 
of the general class of solutions (|l~5|). It is a combination of both monopole and dipole terms 
with the identification 2M = -Mi = -4A 2 M i, so that the constraint M + Mi/2 = 
which guarantees the regularity of the axis is automatically satisfied. This is consistent with the 
interpretation of the solution ( |l4|) given in JlJ]. Indeed, the metric functions p and A given by 
di~4| ) can be written as /i = /^o + A*i and A = Ao + Ai + Aoi, where /io, Ao are given by di~0|), 
//i, Ai are given by (ph, and A i = M 2 ia(l - a 2 )/4A 2 R 3 . 



4 Limiting cases of the above metrics 

Cases 3a and 2a 

A particular limit of case 3a above in which mi = and A\ —* oo was investigated by 
Bicak and Schmidt || in 1989. In this limit, there are only two particles of vanishing mass and 
their accelerations become infinite. In addition, it is necessary that the parameter m = m\ is 
scaled to zero in such a way that the "monopole moment" Mq = —4mA remains constant. The 
resulting space-time is flat everywhere except on an expanding sphere and on the strut between 
the particles. It therefore describes an expanding spherical impulsive gravitational wave that is 
generated by the two particles which move apart at the speed of light in a Minkowski background 
and are connected to each other by an expanding strut. 

An analogous situation in which (for r > 0) two null particles recede from the origin with the 
speed of light and are connected to infinity by semi-infinite strings was constructed explicitly 



by Gleiser and Pullin 15] in 1989. This describes the space-time representing an expanding 
impulsive spherical gravitational wave generated by a snapping cosmic string. Such a solution 
can also be obtained as a limit of case 2a above |l(|. Again, it is necessary that m is scaled to zero 
such that Mo remains constant. However, as pointed out by Bicak |lq| , the complete solution 
rather describes two semi-infinite strings approaching at the speed of light and separating again 
at the instant at which they collide. 

In the remainder of this section, we will investigate similar null limits of all the remaining 
cases described in sections 2 and 3. 

Multipole particles 

We may first consider the null limit of the general class of solutions (jll-16). These represent 



two accelerated particles with arbitrary multipole structure possibly with strings between the 
particles or with strings connecting each particle to infinity. Here, we let A —* oo while all the 
multipole moments M n are kept constant. In this limit we obtain 
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In this limit, \x in (15) approaches the constant C. Also the first term for A vanishes, and only 
the contribution I = in the second term remains finite. Consequently, the metric functions for 
the solutions which we obtain in the above null limit can be written as 



M = C, 

A = ^-si g n(p 2 + C 2 -r 2 )f: (21) 

n=0 

For the particular values of the constants C and D given by (|l~7|) , this solution appears 
to be a generalisation of the null limit of case 2a which describes a situation of a snapping 
cosmic string. The ends of two semi-infinite cosmic string move in opposite directions with the 
speed of light, generating an impulsive spherical gravitational wave. However, in this limit, the 
multipole structure of the initial particles disappears and the solution is characterised by the 
single constant M = J2ri=o M n /2 n . Thus, the null limit for any accelerating particle with a 
multipole structure is identical to that for a particle with just a monopole term Mq = M. 

Considering the alternative values of the constants given by (^), we obtain a null limit in 
which a spherical impulsive wave is generated by a finite string whose length is expanding at 
the speed of light. However, at the ends of the strings the arbitrary multipole structure of the 
initial "particles" described by the multipole moments M n again disappears. This null limit is 
thus identical to that of the case 3a above as originally obtained in ||]. 

Performing the well-known transformation (see e.g. 

p=^(v — u), t = ±\{v + u) coshx, C = \{ v + u ) sinhx, (22) 

we can put the solution ( |2l| ) into the standard form of boost-rotational symmetric metrics with 
null coordinates u and v, 

ds 2 = e A dudv - \{v - u) 2 e - " dcf) 2 - \{v + u) 2 e^ d X 2 . (23) 

Since p 2 + C, 2 — t 2 = —uv, the metric functions take the following form: 

00 M °° M 

snapping string : /j, = - ^ — ^ , A = [Q(uv) - Q(-uv) } ^ — ^ , 

n=0 ra=0 

00 M 

expanding string : fx = , A = 2 @{uv) ^ — — , 

n=o 2n 

where is the Heaviside step function. Notice that in both cases /x is a constant, but there 
is a discontinuity in the otherwise constant value of A with the step 2 Y^=o M n /2 n on the 
null cone uv = 0. However, it is possible to find a transformation to coordinates in which the 
metric is continuous everywhere. In the region where the functions \x and A are constant, the 
transformation 

U = \ue x+ »l 2 , V = ±ve-^ /2 , ^ = Xe M , (24) 
brings the line element ( f23|) into the form 



ds 2 = 4dl(dV-A 2 d(j) 2 -B 2 dip 2 , (25) 

where 

A = V -Ue~ ix+ ^ , B = V + Ue- {x+ ri . (26) 

This metric is explicitly continuous, including on the null cone given by U = 0, even if there 
has been a discontinuity in A at uv = in the original coordinates. Such a step is removed by 
the compensating discontinuity in the expression for U in (E 



The solution for a snapping string can thus be written in the continuous form (|25|). By 
substituting the corresponding A and \i into (p6|), and restricting to the region v > 0, we obtain 
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(27) 



The metric (pq) , fl27| ) is exactly of the type constructed previously by Gleiser and Pullin [15| by 
a different method. It represents an impulsive spherical gravitational wave propagating in the 
Minkowski universe. However, outside the wave (where IA < 0), there extends the "snapped" cos- 
mic string which is characterized by a deficit angle (1 — /3)2tt, where here (3 = exp(^^ M n /2 n ). 
This generalizes the previous result presented in which was obtained as the null limit of case 
2a describing two accelerated monopole particles. And, as in this particular case, the complete 
solution has to be extended symmetrically to negative times, as has been arg ued in @]. 
Similarly, we can write the solution for an expanding string in the form (p5|): 
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9(W)exp -2£ -p + e(-w) 
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V Mr, 



6(W)exp -2^ +e(-w) 
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(28) 



In this case there is a string with the deficit angle (1 — f3~ 1 )2n in the Minkowski space inside 
the impulsive wave in the region U > 0. 

There exists an alternative continuous form of the above solutions. By applying the trans- 
formation ([24|) followed by U = —2U, V = V, Z = (ip + i <fi) , we may convert the metric 
(|23| ) with constant A and fi, in the region u < 0, v > 0, into the form 



ds z 



-2 I VdZ + UHdZ \ 2 -2dUdV , 
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where 77 = — ie and {7 > 0. Performing now a different transformation 
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of the metric (|2"3|) in the region u > 0, v > 0, where J7 < 0, we obtain 



ds< 



-2V 2 dZdZ - 2dUdV . 



(31) 



It is now obvious that we can match the two metrics ( p9[ ) and (31) across the null cone u = 0, 
which corresponds to £7 = 0, so that the resulting metric 



ds 2 



-2\VdZ + Ue(U)HdZ\ 2 -2dUdV , 



(32) 



is continuous. This metric is contained within a general class of metrics which describe impulsive 
spherical gravitational waves (see §, 0, @ and |]). In the general case, the function H(Z) 



can be obtained as the Schwarzian derivative of a "warp" function h(Z) which permits a geomet- 
rical interpretation |19[| . The above solution for a snapping string is given by h = e^ z , where 
P = y/2 exp(££° =0 M n /2 n ) corresponding to 7f = -\ exp(2 ££° =0 M^/2 71 ). The alternative case 
for an expanding string is given by P — ► corresponding to 77 = — | exp(— 2^^L M n /2 n ). 



An alternative limit and Case 4a 



It can immediately be observed that if the condition J2^=o M n /2 n = is satisfied (i.e. when 
the original metric is regular everywhere on the axis except at the locations of the particles), 
the corresponding null limit results in a trivial Minkowski spacetime without an impulsive wave. 
However, it is also possible in this case to perform a more involved null limit which results in a 
different class of spacetimes than those described by (|2l|). The sum in the second term for A in 
( |l5| ) can be written as 

ElrE^y^Elr.^E^.E^y.El?. (33, 

n=0 1=0 x y/ n=0 v y/ n=\ 1=2 x y/ n=l 

When X^nLo -^n/2 n = 0, the second term in the expansion 

2a \^Mn = 2M a ^ 



x — y J , 2 n x — y 



becomes dominant, and a different null limit can be obtained. Considering that x - y ~ A 2 m 
the null limit A — > oo, we may here rescale the remaining monopole moment M$ in such a way 
that M = MqA~ 2 remains constant. With these assumptions and using ( p0| ) we can write the 
final form of the non-trivial null limit as 
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P 2 + ( 2 

9 2 (C 2 -T 

(p 2 + C 2 -r 2 ) 4 ' (p 2 + C 2 -r 2 ) 2 



A = —M 2 , P 2{C \ - A +M - f P l ■ 5' + T ' sign ( p 2 + C 2 - t 2 ) . (35) 



This represents a unique alternative solution which can be obtained in the null limit described 
above if J2n=o M n /2 n = 0, after a suitable rescaling of the multipole moments. However, even 
in this more involved limit, the multipole structure is again completely "erased". Moreover, the 
metric is not acceptable as a description of an impulsive gravitational wave since it is singular 
on the null cone p 2 + Q 2 = t 2 . It may further be noted that when Mq = 0, the solution is 
simply Minkowski space. 

Notice finally, that the solution (|35|) exactly corresponds to that which is obtained in the 
null limit A — > oo from the metric QliD of case 4a with the identification M = —2Mq\. 

Cases 2b, 1, 2 and 3 

Concerning case 2b, it may immediately be seen from ( |l3| ) that the null limit in which 
A2 — > 00 does not exist. The terms in the metric functions which represent the external field 
responsible for the acceleration diverge. 

We have now considered all the above special cases which involve just two accelerating 
particles. We therefore return to the remaining cases which involve two pairs of distinct particles. 
In all these cases, it is only possible to consider null limits as A% — > 00 (since A2 > A\). In cases 
1 and 2, this limit diverges (p ~ A 2 .). However, in case 3, the limit for p is given by 

2mi 1 

fi = 



AtRx A\\p 2 + C 2 - r 2 |' 

and this can be seen to diverge on the expanding spherical surface p 2 + £ 2 = t 2 . Null limits of 
the cases 1, 2 and 3 must therefore be considered to be unphysical. 



Case 4 



It now only remains to consider the limit of case 4. We again consider the limit as A2 — > 00, 
but we can now do this while scaling m,2 to zero in such a way that A2rri2 = —A\mi = constant. 
In this case the limit is given by 

2m 

M = 



A 



AR 

2 

m 



t 2 ) - sign(/> 2 + C 2 - t 2 ) 4mA, 



(36) 



A 2 R ^ v* ■ ' ' R 

in which R = R\ is given by (jllj) , and the subscript 1 has been omitted from mi, A\ and R\. 
This solution contains two free parameters m and A (although it is not possible to consider the 
limit as A — > while keeping m finite). It represents a snapping string of finite length whose 
outer ends are accelerating apart as shown in figure 10. 

T 




Figure 10. A space-time diagram for a finite snapping string 



The two outer particles of mass m are caused to accelerate by the strut or string (there is 
a deficit angle if m < 0) connecting them. However, in this case, the finite strut breaks (or the 
string snaps) at its midpoint and the two broken ends separate at the speed of light. 

Apart from the step change in A, the metric ([36 ) is identical to that of cases 2a and 3a 
as given in (|lC|). For the general family of metrics (l|), the curvature tensor components are 
linear in A, and its discontinuities give rise to impulsive components on the null cone. In fact, 
the Ricci tensor components still vanish everywhere. However, an impulsive component arises 
in the Weyl tensor on the null coneQ. This is interpreted as describing an expanding spherical 
impulsive gravitational wave that is generated by the snapping of the string. 

The complete solution of course is time symmetric as illustrated in figure 10. This situation 



resembles that of the known solutions for snapping infinite cosmic strings [15[, [16]. However, 
this solution differs in that it is not flat on both sides of the impulsive wave. 



5 Conclusion 

We have reviewed the Bonnor-Swaninarayan family of boost-rotationally symmetric solutions 
which was generalised by Bicak, Hoenselaers and Schmidt. We have written these solutions in 
a unified way paying particular attention to certain limiting cases. These solutions describe the 
accelerated motion of pairs of particles possibly attached to conical singularities on the axis of 
symmetry. For the case of multipole particles, we have presented a new simpler form of the 
solution. 

We have specifically investigated the possible null limits of these solutions. In most cases, 
the limits are trivial or not physically acceptable. However, in some cases we have obtained 
physically interesting limits which describe snapping or expanding cosmic strings generating 
spherical impulsive gravitational waves. For multipole particles, we have shown that the internal 
multipole structure vanishes in this limit, leaving the same solution for a snapping or expanding 
cosmic string. We have also presented a new solution for a snapping cosmic string of finite 
length. 

1 These statements apply also to the null multipole limits (EHl) as well as to the null limits of cases 2a and 3a. 
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Appendix 

The solution ([l5|) for multipole particles was first presented in [14] using prolate spheroidal coor- 
dinates. The expression for /i and the first term for A in ( |i"5| ) are known from the standard Weyl 
solutions. However, the second term in A arises as a 'mixed' term representing the interaction 
between the multipole terms and the 'boost potential'. Denoting the interaction term for the 
n th multipole by X n , this is required to satisfy the equation 

2(n + l) /Pn+l _ 2P n 
v 7 /? 2 - 4a/3 + 4 \f3 n+1 p 



d P X n ~ Jm \ n , =T ( TmTT ~ 1^+2 ) > ( 37 ) 



(this is a correction of equation (3.11) in fl4|| in which ^A n — 7 mjn — 

p-(n+l)p n j where the 

independent coordinates are a = (1 — xy)/(x — y) and (3 = x — y, and the Legendre polynomials 
have argument a. We have found that this can be integrated directly by introducing the function 

A n = A n — ^A n _i. (38) 



Using the recurrence relation for Legendre polynomials and (37), we find that A n must satisfy 



which can immediately be integrated. Setting the integration constant to zero and using a 
recurrence formula derived from ( |38[ ) with Ao = —fi~ X \J ft 2 — 4a/3 + 4, we finally obtain the 
integral 

1 / n /2\ l 

^ = - ¥ - 4a/3 + 4 g i^-j P l . 

In prolate spheroidal coordinates, this is 

1 /x + y x n / o \/ 



2" - y/ ^ - yj 



which, with the coefficient M n , gives the second term in (fi5|) above. 
The equivalent expression presented in 14], equation (3.13), is 



2" \x-y 

where 



(n + l)(P n / n+1 -P n+1 / n ), (40) 



fe=0 1=0 



Note that in equation (3.13) of 14], there is also a term P n / n+1 (oo) — P n+1 / n (oo), where 
/ n+ i(oo) = J2k=o kTi Pn - kPk - However > usm S tne identity (n + l)(P n+1 Q n - P n Q n +i) = 1 
and the standard definition Q n (a) = ^P n (a) ln[(l + a)/(l — a)] — f n (oo), it can easily be 
shown that (ra + l)[P n / n +i(oo) — P n+ i/ n (oo)] = 1. This term thus represents only an additive 
constant which can be omitted. 

Since both the solutions (|3^) and ( f40"P satisfy the same differential equations, they can differ 
only by a constant. However, by considering the particular value at x = 1, y = 0, it can be seen 
that this constant must be zero. From this we may now deduce the relation 

n ( 2 \ l 
(n + l)(PJ n+ i - P n+l f n ) = £ Pi . 

By expanding the left-hand side in powers of 2/{x — y) and comparing coefficients, we obtain 
the following non-trivial identity for the Legendre polynomials. 

Theorem: For non-negative integers m and n, with m < n, define the functions X^ 1 (a) by 

n 

x {m) _ _J_ p p 

n — / j k+1 n—k 1 k—m > 
k=m 

where P n (a) are Legendre polynomials, and X^\ = 0. Then 

P n Xn ^ — P n +iX^ l _ 1 = ^j-j- P m (41) 
for arbitrary m and n and for arbitrary value of the argument a. 

The left-hand side of the above identity is a polynomial of order 2n — m. It can therefore be 
expanded as a series of Legendre polynomials. However, remarkably, all the coefficients vanish 
identically except for one. 



